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An Abelian two-Higgs model of strongly correlated electrons: 
phase structure, strengthening of phase transition and QCD at finite density 
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We investigate non-perturbative features of a three-dimensional Abelian Higgs model with singly- 
and doubly-charged scalar fields coupled to a single compact Abelian gauge field. The model is 
pretending to describe various planar systems of strongly correlated electrons such as high- Tc super- 
conductivity in the overdoped regime and exotic materials possessing excitations with fractionalized 
quantum numbers. The complicated phase structure of the model is studied thoroughly using nu- 
merical tools and analytical arguments. In the three-dimensional space of coupling parameters we 
identify the Fermi liquid, the spin gap, the superconductor and the strange metallic phases. The be- 
havior of three kinds of topological defects - holon and spinon vortices and monopoles - is explored 
in various phases. We also observe a new effect, the strong enhancement of the phase transition 
strength reflected in a lower order of the transition: at sufficiently strong gauge coupling the two 
second order phase transitions - corresponding to spinon-pair and holon condensation lines - join 
partially in the phase diagram and become a first order phase transition in that region. The last ob- 
servation may have an analogue in Quantum Chromodynamics at non-zero temperature and finite 
baryon density. We argue that at sufficiently large baryon density the finite-temperature transi- 
tion between the (3-flavor paired) color superconducting phase and the quark-gluon plasma phases 
should be much stronger compared with the transition between 2-flavor paired and 3-flavor paired 
superconducting phases. 
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I. INTRODUCTION 

Gauge models involving multiple scalar fields coupled 
to an Abelian gauge field are applicable to a large va- 
riety of systems such as multi-band superconductors^, 
liquid metallic hydrogen^, easy-plane quantum antifer- 
romagnets^, etc. These models have interesting phase 
structure and are distinguished by a copious zoo of topo- 
logical defects. Usually, all scalar (Higgs) fields in these 
models are considered to be alike so that they are all 
equally charged and minimally coupled to the gauge field. 
Contrary, in this paper we consider the two-Higgs model 
with a gauge field and two Higgs fields with unequal 
charges. Our study is motivated by the fact that the 
charge-asymmetric two-Higgs gauge model can emerge 
as an effective description of unconventional supercon- 
ductivity^. 

Despite the two-Higgs model is formulated in a very 
simple way, it can actually capture basic properties of 
various systems: instanton (monopole) plasma described 
by the compact Abelian gauge model, superfluidity (the 
XY spin model), the Z2-gauge model which has ana- 
logues in particle physics and the coupled Z2-gauge- 
XY model describing fractionally charged excitations in 
strongly correlated electron systems. It shares also a sim- 
ilarity with a Ginzburg-Landau model with two vector 
fields^ which was suggested to describe extended s- and 
d-wave superconducting granular systems. 

The common feature of all unconventional supercon- 
ductors^ with high critical temperatures is the presence 
of copper oxide layers. The layered structure is seen 



both in electronio^ and optical^ anisotropic structure of 
the cuprates. The anisotropy is a key ingredient of var- 
ious approaches to this phenomenon^. The Cu02 layers 
are associated with conductance plates while the atoms 
in the space between the layers form a so-called charge 
reservoir which supplies the charge carriers to the planes. 
The charge reservoirs themselves are almost insulating in 
the superconducting phase^. The charge carriers can be 
either electrons or holes depending on the nature of the 
dopant. The fraction of the carriers p in the planes is con- 
trolled by the doping level x which is usually encoded in 
the chemical formula of the cuprate oxides (i.e., p = x in 
the structurally simple La2_KSr2; CUO4 as a prototype of 
many of the cuprate materials). In the clean limit, x — 0, 
the cuprates are Mott insulators, while at certain x the 
cuprate becomes a poor conductor which - at low enough 
temperature - turns into a superconductor. Nowadays, 
the critical temperatures have climbed the level of 140 K 
in Hg-based cuprates. 

We concentrate on the in-plane mechanism of the high- 
Tc superconductivity restricting ourselves to the slave- 
boson approach in the t-J model. That model is used 
to describe the ground state of the high- Tc superconduc- 
|-Qp9 ,io,ii charge carriers (electrons or holes) in the two- 
dimensional copper-oxide plane. The t-J Hamiltonian^ 
describes hopping holes (or electrons) and localized spins 
in a plane: 
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n.rij) . (1) 
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The first term specifies holes (electrons) moving with- 
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out flipping the spin cr. Double occupancy is explicitly 
forbidden by the presence of the projectors 

Pij.cr = (1 - - n-j,^) . 

The second term describes the anti-ferromagnetic Heisen- 
berg coupling between spins located at the copper sites. 
Here 

is the spin operator, cj^, Cia are the hole (or electron) 
creation and annihilation operators, and 

denotes the occupation numbers. 

According to a popular scenarioi^, the electron degree 
of freedom can be split into the spin and charge con- 
stituents (spinon and holon, respectively). The splitting 
gives rise to an internal gauge degree of freedom, with re- 
spect to which the spinous and holons have positive and 
negative charges (say, and —1, respectively). The 
internal group is necessarily compact and this leads to 
a specific interaction between the spin-charge separated 
constituents to be discussed later. 

Under certain conditions the spinon particles get 
paired and form pairs similarly to Cooper pairs in or- 
dinary superconductivity. Then the pairs of spinons are 
presented by a spinon-pair field which has charge +2 with 
respect to the internal group. Therefore, in a mean field 
approach, the system is described by two scalar (Higgs) 
fields: the holon and the spinon-pair field with internal 
charges —1 and +2, respectively. Both kind of fields in- 
teract via the exchange of an internal gauge field which 
is compact by construction. 

It is important to note that the group for the inter- 
nal gauge degree of freedom has nothing to do with the 
usual Maxwell electromagnetic group. For example, the 
internal degree is compact while the electromagnetism is 
described by a non-compact Abelian group. The original 
spinon is an electromagnetically neutral excitation while 
the holon is the only constituent which carries the elec- 
tric charge. Thus the charge carrier is the holon while 
the spinon may affect the properties of the strongly cor- 
related material only indirectly: the formed spinon-pairs 
interact via the compact gauge field with the holons. 

Entering a more technical description, the creation op- 
erators are decomposed a a^°i^^ 

4 = flb^ , (2) 

with the constraint 

4/n + 4/a + blh = 1 . (3) 

Here fi^ is a spin-particle ( "spinon" ) operator and bi a 
charge-particle ("holon") operator. 



In addition to the ordinary electromagnetic (external) 
gauge symmetry, 

^(f)cxt ■ Cjfj > e , fia ^ fi<77 ^ ^ ^icr (4) 

the spin-charge separation naturally introduces an (in- 
ternal) compact U{1) gauge freedom, 

^(f)int • ^i<7 ^ Cicr , fifj > € fia: ^ ^ ^icr (5) 

which plays an essential rol o^^i^^ in understanding the 
physics of strongly correlated electrons. The spin-charge 
separation idea may also be applied to various systems 
including the general case of non-relativistic electrons — 
as well as the case of strongly interacting gluons in Quan- 
tum Chromodynamicsi^. 

The effective theory of superconductivity can further 
be simplified and reformulated in terms of lattice gauge 
models^ii^ii^ii^iii^, see recent reviewsi^. Thus, the t-J 
model ^ is related to a compact Abelian gauge model 
with the internal symmetry ([5]), which couples holons and 
spinons. As in usual BCS superconductivity, under ap- 
propriate conditions the spinons couple and form bosonic 
quasiparticles. In a mean field theory one can define fields 
which behave under the gauge transformations ([5]) as: 

- (/.T/,i-/a/.-T)-A,,.e^("'+"^). (7) 

The phase of the field x represents nothing but the com- 
pact U{1) gauge field, 

% = argx^i ^ % + (da),^- 

with (da)j^- — aj —ai, and the radial part, x = KXy)!; is 
the so-called "resonating valence bond" (RVB) coupling. 
The doubly-charged spinon-pair field A is analogous to 
the Cooper pair. 

At high temperature the RVB coupling vanishes, x = 
0, and the system is in the Mott insulator (or "poor 
metallic" ) phase. With decreasing temperature x ac- 
quires a non-zero value, eventually enabling the forma- 
tion of a spinon-pair condensate A = |(Aij)| and/or 
of a holon condensate b = (bi)^. Therefore, four 
phasesi^iiiii^ may emerge: 

(i) the Fermi liquid (FL) phase with & 7^ 0, A = 0; 

(ii) the spin gap (SG) phase with 6 = 0, A 7^ 0; 

(iii) the superconductor (SC) phase with 6 7^ 0, A 7^ 0; 

(iv) the strange metallic (SM) phase with 6 = 0, A = 0. 

The ground state of the superconducting layer is pro- 
posedi^iii to be described by a compact Abelian two- 
Higgs model (cA2HM) in three dimensions with a [/(I) 
gauge link field 9, a singly-charged {qh = 1) holon field 
and a doubly-charged {qs = 2) spinon-pair field ^s- 
More precisely, the classical three-dimensional statistical 
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Insulator 
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Superconductor 



Fermi liquid 

(b^O, A =0) 



Dopant concentration, x 

FIG. 1: The phase diagram in the plane "temperature-dopant 
concentration" (suggested in Ref . Ha) . 



model describes the ground state of the zero temperature 
fwo-dimensional spin-charge separated quantum system 
of holons and spinous coupled together by the internal 
compact gauge field. The model can be considered as 
a phenomenological extension of the Ginzburg-Landau 
model for ordinary low- Tc superconductors. 

The model is applicable to the overdoped regime of 
high- Tc materials (called sometimes the SM region^'^^) 
where the SU (2) particle-hole symmetry is explicitly bro- 
ken. Further arguments for justifying this approach and 
a discussion of its limitations can be found in Ref. [J. 

Since high- Tc materials are type-2 superconductors, we 
restrict ourselves to the London limit in which the ra- 
dial parts of both Higgs fields = |$fe| e^'^'', k ^ h,s 
are frozen, \^h,s\ = const. The action of this compact 
Abelian two-Higgs model is 



5'cA2HM = X! 



-Kh ^ cos{diph + d)i - ^ cos{Aips + 2e)i , (8) 
; / 

where Op = [dO) p is the standard lattice plaquette. 

The model (|8]) obeys a lattice version of the U{\) in- 
ternal gauge symmetry (O: 



iPk ~ qkOi : 



(9) 



with k = h,s. It describes the hole (electron) "con- 
stituents" by the dynamical holon ipf^ and spinon ip^ 
phases, which strongly interact via the dynamical gauge 
field 6. The inverse gauge coupling /3 in Eq. ([8]), /3 = 
X/ + Xbj is given by the diamagnetic susceptibilities Xf 
of the spinon and Xb of the holon fields. The hopping pa- 
rameters Kfc are connected to the doping concentration x 
and the couplings t and J of ([1]) as follows: Kh (x t ■ x 
and Kg (X J (seek^. 

The phase diagram of the model ^ and the basic prop- 
erties of the topological defects were studied for a limited 
set of coupling parameters in our preliminary investiga- 
tion^^. The aim of the present paper is to extend the 
study to a much larger coupling range using both nu- 
merical and analytical tools. We also identify the order 



of various phase transitions and confirm the existence 
of a phase transition strengthening effect which was first 
suggested in Ref. [2l|. The enhancement is reflected in de- 
creasing the order of the phase transition at a common 
joint segment of two second order transitions. The sig- 
natures of unexpected strengthening of the phase transi- 
tion were later observed in a different model describing a 
gauge field coupled to two Higgs fields of equal charges^^. 

The structure of the paper is as follows. In Section |TT] 
we show that the model contains three types of topolog- 
ical defects: two types of vortices and one type of mono- 
poles. We discuss the simplest features of the topological 
defects and derive their effective action using analytical 
tools only. In Section Hill we describe the phase diagram 
of the model in the three-dimensional space of gauge, 
holon and spinon couplings using known results avail- 
able for less complicated systems. The limiting cases of 
the three-dimensional phase-cube are analyzed in detail 
and the possible structure of the phase transition sur- 
faces in the SD-coupling space is pointed out. Section HVl 
is devoted to a numerical investigation of the phase di- 
agram by Monte-Carlo simulations. We analyze various 
two-dimensional cross-sections of the 3D-phase diagram 
identifying numerically the order of the phase transitions. 
The behavior of thermodynamical quantities in combina- 
tion with that of the densities of the topological defects 
allows us to identify the nature of the phases in different 
regions of the 3D-coupling space. In the same Section we 
confirm the effect of the phase transition strengthening 
due to merging transition lines. The phase transition en- 
hancement may also be relevant for Quantum Chromody- 
namics which describes the theory of strong interactions. 
We point out in Section |V] that the phase diagrams of 
the cA2HM and QCD contain common features includ- 
ing the joining of transition lines. We suggest that QCD 
at high temperature and high baryon density may expe- 
rience the same strengthening effect. The last Section is 
devoted to our conclusions. 



II. PROPERTIES OF TOPOLOGICAL DEFECTS 

The model ^ contains three kinds of topological de- 
fects: a monopole and two types of vortices, referred to 
as the holon and the spinon vortejs^ii^. The monopole 
has magnetic charge 27r?l while the holon (spinon) vor- 
tex carries magnetic flux quanta 2TTh (vrft) of the gauge 
fleld 9. One monopole is simultaneously a source of one 
holon vortex and two spinon vortices. In this Section we 
discuss some basic properties of these defects. 

Despite the model ([8|) is formulated in the Wilson 
representation (with an action of cosine- type), the ba- 
sic properties of the topological defects can be guessed 
from the so-called Villain representation^^ which is more 
suitable for analytical considerations. A similar set of 
transformations was performed in a non-compact model 
in Ref. L. The principal difference between the compact 
(considered here) and non-compact modelsi is the pres- 
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ence of monopoles, and, as a consequence, a richer phase 
diagram due to the existence of confining phases. 

The ViUain representation of the partition function of 
the model ^ is 

^ = / E / D^h / D^, E E 

exp{-/3||d0 + 2^n||' - K,.||d^^ + 9 + 2TTl,,\f 

-iis\\d^s + 20 + 2TTls\fY (10) 

where the ViUain coupUngs and kg correspond to the 
Wilson couplings and k^. There are three integer- 
valued forms: the plaquette form n and two link forms 
Ih and Is. 

The definition (|10p is written in a convenient con- 
densed form by means of differential forms on the lat- 
tice^. In brief, the notations are as follows. Let a and b 
be two r-forms on the lattice. Here r = corresponds to 
scalars (with the support on sites, cq), r = 1 corresponds 
to vectors (with the support on links, ci), etc. Then the 
scalar product (a, 6) is defined as a scalar product over 
the whole support of the r-form (sites, links etc.) over 
the lattice. Thus, for two vector forms ("one- forms"), 
we have (a, 5) = ^jajfoj. The modulus squared is then 
defined as ||a|p = (a, a). The finite-difference operator 
"d" increases the rank of the form by one, r — > r -I- 1 
(thus having the meaning of a gradient), while the oper- 
ator S = *d* (to be used below) decreases the rank of the 
form, r —^ r — 1 (having the meaning of a divergence). 
The duality operator - which switches forms between the 
original and the dual lattices (by element-wise equating 
the values assigned to the dual to each other supports) - 
is denoted as "*" . 

The Berezinskii-Kosterlitz-Thouless (BKT) transfor- 
mation?^ allows us to rewrite the partition function (jlOp 
in terms of monopoles and vortices. The monopoles ap- 
pear due to the compactness of the gauge fields 6, and two 
types of vortices arise from the presence of two indepen- 
dent species of Higgs fields. The monopole "trajectory" 
(in 3D actually located on cubes) is denoted as j, the 
vortex "trajectory" (in 31? located on plaquettes) are ah 
(the holon vortex) and cr^ (the spinon vortex). Using the 
standard approach we represent the integer-valued form 
n as a sum n — dq + m[j] of the co-closed surface dq and 
the surface m[j] spanned on the monopole trajectory j: 
*j = S*m[j], or, equivalently, j = dm[j]. The two-form 
n can also be represented in Hodge-de-Rham form as a 
sum of a closed and a co-closed part 

n ^ 5A-^j + d{A-''^ 5m[j] + s) , (11) 

where A = d(5 + (5d is the lattice Laplace operator, A^^ 
is its inverse, and s G Z(ci), m S ^(02), and j S ^(03) 
are integer-valued forms. The monopole current is closed 

6*j = . (12) 



Substituting Eq. (Ill|) into the first term under the ex- 
ponent in Eq. ^TU\i we get 

de + 2Tm^ dA + 2T:5A-^j , (13) 

where 

A = e + 2T:{A-^Sm[j] + s) (14) 

is a non-compact gauge field. 

The same trick can be performed with the phases iph 
and ips (in the second and the third term in the exponent 
in Eq. (Hn}) 

Ik = -^A-V'"^ + d(A-Mpfe[f7f ] + r,) , (15) 

where cr^"^ G Z(c2) is the closed part of the holon (if 
k = h) ov spinon (if k ~ s) vortex trajectory, 

(5V^°^ e Z(c2) = 0, k = h,s, (16) 

Pk S Z(c2), and Tfe G Z(co). The analogue of Eq. (fT^ 
can be written in the following form 

d<y9fc +qkO + 2nlk = dtpl ''- + qkA + 27r(5A~ . (17) 

Here the vortex trajectories are given by the combina- 
tions 

f^fc = ^k'^ - Ik mij] , (18) 
and the non-compact Higgs phases are 

^r-='Pfe + 27r(A-i<5pfe[ar]+r,). (19) 

The vortex trajectories end on the monopole "trajecto- 
ries" , which are points (instantons) in three-dimensional 
space-time, 

S*(Tk+qk*.j = 0; k = h,s. (20) 

The above equation means that there are one holon vor- 
tex and two spinon vortices - each carrying a correspond- 
ing elementary flux - attached simultaneously to each 
monopole. A typical configuration monopole-vortex con- 
figuration is shown in Figure [21 



Spinon-pair vortex 




Spinon-pair vortex 



FIG. 2: An example of an monopole-vortex configuration. 

Substituting Eqs. ^ and ^ into ^ and perform- 
ing the Gaussian integration over the non-compact fields 
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A, (fi'^'^', and (p'^''^', one gets the BKT representation of 
the two-Higgs model: 



E 



(21) 



<5*o-fc+9fc *J=0 

where an irrelevant constant factor is omitted. The mo- 
nopole and the vortex actions are, respectively, 

1 



An- 



An- 



A- 



E 

k.k'—h.s 



Cfe, ICkk' CTk' 



with the differential operator /C 



A 



'fcfc' 



qk> A + n-ii 



(22) 
(23) 

(24) 



In Eqs. (p2ll . ([24| denotes the mass of the photon, 



- 4k,) , 



and oJk are the normalized weights 



(25) 



(26) 



The (self-)interaction of the monopoles and vortices 
can readily be read off from Eqs. ^T^, and (^5)) . 
The monopoles interact via a massive photon exchange, 
and thus the interactions of the monopole "trajectories" 
are suppressed at large distances. 

The vortex interactions contain a long-range term. 
The origin of the long-range forces is simple: in the 
Abelian model with one Higgs field the massless Gold- 
stone mode is eaten up by the longitudinal component 
of the gauge field. Thus, the gauge field becomes mas- 
sive, whereas the massless Goldstone boson disappears. 
This is not the case in the present model with two Higgs 
fields: one of the Goldstone bosons can be absorbed into 
the longitudinal component of the gauge field, while the 
other remains alive. Thus, the two-Higgs system with 
one gauge field will always have one massless excitation 
which, in particular, leads to a long-range interaction be- 
tween the vortex "surfaces" . 

The vortex- vortex interaction term (|24p can also 
be rewritten in the form 



■s: 



ST[a,,a,] (27) 



where the holon-holon, spinon-spinon and holon-spinon 
vortex interactions are, respectively. 



■'hh 



An Kh (Th 



1 



A 



A- 



1 



1 

A 

1 - 
A. 



A-|-to2 



1 

A 



(28) 
(29) 
.(30) 



From these expressions it can be seen that the short- 
range interaction between the parallel segments of all 
vortices is always repulsive. However, the presence of 
the long-range component acts in a different way: seg- 
ments of equal type vortices are repulsive (holon-holon 
and spinon-spinon), while the parallel vortex segments 
of different types (holon-spinon) are always attractive. 
The respective strengths of the repulsion and attraction 
depend on the weights of the holon and spinon vortices, 
Eq. (pS)) . which, in turn, depend on the strengths of the 
holon and spinon-pair condensates, Rk oc 

The presence of the massless mode in the interaction 
between vortices (^5)) . (f^. ([50)1 on the Lagrangian level 
does not mean that on the quantum level the interac- 
tions remain unscreened. In contrary, we expect that the 
massless mode should disappear due to non-perturbative 
effects. This happens, for example, in the monopole gas 
of the three-dimensional compact Abelian model^^^. In- 
deed, in this example the bare interaction between the 
monopoles is of the Coulomb type, (j, A~^j), while all 
correlations in the statistical ensembles of the monopo- 
les are exponentially suppressed at large distances by a 
Debye mass which appears due to monopole interactions 
in the plasma regime. A similar effect is expected in the 
considered system of Coulomb-interacting vortices. In 
the dilute ensemble of vortices the screening mass is ex- 
pected to be small compared to the mass of the photon 
TO, Eq. (|25p . but still non-zero. Therefore one may ex- 
pect that on the quantum level the Goldstone mode may 
disappear. 



III. THE PHASE DIAGRAM 

The internal structure of the three-dimensional phase 
diagram of the (London limit version of the) compact 
Abelian two-Higgs model dS]) is rather complicated. How- 
ever, the faces and the edges of the cube representing the 
"compactified" phase diagrams can be drawn relatively 
easily because they are related to various well-known 
condensed-matter systems. These limiting cases of the 
phase cube correspond to appropriate combinations of 
vanishing and/or infinitely large couplings /3, and k,. 
Below we discuss these limits in detail. 



A. The (/3, Kh)-faces 

The Kg parameter defines the coupling of the spinon- 
pair field to the compact gauge field. Below we con- 
sider two limiting cases corresponding to vanishing and 
infinitely strong coupling k,. 



1. The K,s = G face: Q = 1 compact Abelian Higgs model 

The Kg — face of the cA2HM corresponds to the Q = 
1 Higgs model with a compact Abelian field (cAHMq=i) 
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in three dimensions: 



3cAHM 



Q=i 



— — /3 cos ( 



K/l ^ C0s(d(/3ft 

/ 



^);.(31) 



On that face the holon condensate is coupled to the com- 
pact J7(l) gauge field while the spinon-pair field is decou- 
pled from all other fields. The phase of the spinon-pair 
field is disordered which implies condensation of the spi- 
non vortices and, therefore, vanishing of the spinon-pair 
condensate A. 

The cAHMq=i has extensively been studied in the 
literature and it seems that a consensus on the phase 
structure is reache d^^'^^i^^'^° . The phase diagram in the 
(/3,K?i)-plane is plotted in Figure [H 




FIG. 3: The phase diagram on the Ks = face corresponds to 
the compact Abelian Q = 1 Higgs model with the action (I31|l . 
The SM and FL phases in the interior of the diagram, as 
well as the limiting models (XY , cQED and the trivial cases) 
are explicitly indicated. The Xy-critical point is given in 
Eq. 



It contains two phases: 

(i) the SM phase at small /3/small kk, 

(ii) the FL phase at large /3/large k^. 

The FL phase is the broken/Higgs phase with nonzero 
holon condensate 6 7^ 0, and the SM phase is the con- 
fining/symmetric phase with 6 « 0. The distinction be- 
tween these phases is blurry since the broken phase is 
always partially confining while (at the opposite end) 
even in the deeply confining phase traces of the Higgs 
condensate can be found. The usual order parameter - 
the Higgs (holon) condensate - is, strictly speaking, non- 
zero in the whole phase diagram, and thus it is usually 
said that these phases are "analytically connected" . As a 
consequence there is no local order parameter in terms of 
the primary fields of the model, which could in principle 
discriminate between the phases. 



Despite the boundary between broken and confining 
"phases" has definitely not the characteristics of a phase 
transition in the thermodynamic sense, these phases may 
be discriminated by the condensation properties of topo- 
logical defects. A non-thermodynamic boundary of this 
type is known in the literature as a Kertesz line^^'^ which 
is defined as a line where the vortices start /cease to con- 
dense. This line - which is not plotted in Figure [3] - 
has been thoroughly studied in the three-dimensional 
Abelian Higgs model in Ref. [H and was suggested to 
appear in models of particle physics as well^^. Note 
that in the present context the word "line" refers to a 
two-dimensional coupling parameter space, while in, say, 
a three-dimensional parameter space the corresponding 
manifold is a surface. 

The edges of the cube in the Ks = 0-plane can be an- 
alyzed following Refs. H^IUlllli At the k/, = edge 
the model is basically a plasma model. In fact, along 
the K/i = Ks = edge the cA2HM reduces to the pure 
compact ?7(f ) gauge theory. 



5, 



cQED 



cos Vp , 



(32) 



which is known to be confining at any value of /3 due 
to the presence of monopoles-^. The monopoles are in- 
teracting as Coulomb particles, thus forming a magnet- 
ically neutral plasma. In three dimensions the monopo- 
les are pointlike, i. e. instanton-like objects. Along the 
K/i = = edge the vacuum of the model possesses a 
mass gap at any non-zero density of monopoles, which is 
realized at any finite [3. The mass gap is given by the 
Debye mass of the monopole plasma. 

Along the edge /? — > 00, Kg =0 the model describes 
a superfluid. Indeed, along this edge the gauge field be- 
comes constrained to the trivial vacuum, d9 = 0. The 
constraint is resolved as = da. Identifying Lph — ^1 
the model becomes the 3D XY spin model 



S 



XY 



cos{dip)i 



(33) 



where we have set = k^^. The XY model is known 
to possess a second order transition'^^ at the critical point 



« 0.45420. 



(34) 



shown by the dot in Figure [31 

The edges /3 = and ^ 00 correspond to trivial 
theories so that they do not posses any phase transitions. 



2. The Ks 



00 face: XY spins coupled to Z2 gauge field 



In the limit Ks — > 00 the coupling of the spinon-pair 
field and the gauge field is tight. Mathematically, this is 
expressed in the form of the constraint 



29 + dips = 27rTO , m G Z , 



(35) 
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which is to be fulfilled at each link. The constraint p5|) 
has the solution 9 = (d(^s/2 + 7rm)2,r, where m £ Z is 
chosen such that 9 £ (— tt, tt]. Substituting this solution 
back into Eq. ([5]), and introducing the Z2 gauge field 
ai = (—1)™', we get that on the face Kg ^ 00 the model 
(III) reduces to a 3D XY-Z2 model with the action 

SxY-Z2 ^ -P^crp - Hh ^ CT; cos{Aip)i . (36) 
p I 

This model describes a XY-type matter field Lp — iph 
'{'sl^ which interacts via the exchange of a Z2 gauge field 
<Ji . The model has a rich phase structure studied numeri- 
cally in Ref.m. The phase diagram ~ plotted in Figure [4] 
- contains three phases: 

(i) the SM phase at small /3/small Kt, 

(ii) the SG phase at large /3/small k^, 

(iii) the SC phase at large /?/large Kh- 




FIG. 4: The phase diagram on the Ks — > cx) face corresponding 
to the model (|36p . The critical points at the edges and the 
phase transition lines in the interior are discussed in the text. 

As one can see from Figure |4] the phase structure of the 
— > 00 face is much richer than that of the Ks — face. 
Let us analyze the edges of this two-dimensional phase 
diagram. At the ~ edge the model is reduced to the 
Z2-gauge (or, the "gauge Ising") modeP^ with the action 

S^i = -/3j2ap. (37) 
p 

This model has a string-like topological object as well 
as vacuum excitations which are sometimes considered 
as prototypes of, respectively, the chromoelectric string 
and glueball excitations in the strongly interacting quark- 
systems in non-Abelian gauge theories'^^. 



The gauge Ising model (|37)) is known to possess a sec- 
ond order phase transition of the Ising-type at the critical 
point 

lim (3c{Kh) = Pf « 0.7613 . . . (38) 

shown by the dot at k/i = edge of Figure^ This transi- 
tion separates the disordered (small (3) and ordered (large 
(3) phases discriminated by the presence or absence of the 
spinon-pair condensate, A = and A 7^ 0, respectively. 
At the same time, at this axis the holon condensate van- 
ishes due to disorder in the XF- variables (to be discussed 
below) . Therefore the small-/3 and large-/3 phase are iden- 
tified with the SM and SG phases, respectively. 

At the (3 ^ 00 edge the Z2 gauge field is suppressed due 
to the constraint ap = 1, so that the model is reduced 
to the XY model (1551) with a coupling k,^^ = Hh- At 
small values of this coupling, Kh < k-^^ , the holon field 
is disordered by vortices and the holon condensate is ab- 
sent. At large values Kh > k^^, the vortices are dilute 
and holon condensation takes place. These regimes are 
separated by the critical coupling which is marked 
by the dot at the (3 ^ 00 edge of the phase diagram. 
Figure ID 

The interior of the phase diagram is also nontrivial. 
The phase transitions, marking the onset of spinon-pair 
(dotted line) and holon (solid line) condensations, are 
departing from the Kh — Q and (3 ^ 00 edges towards 
the center of the phase diagram. There they meet to- 
gether forming a single transition line where the holon 
and spinon-pair condensations occur simultaneously, Fig- 
ure [H This combined line extends down to lower /3, and 
finally it meets with the /3 = edge of the phase dia- 
gram. It is described by a (modified) XY' model with 
the action 

SxY' = - ^ log cosh[K/i cos(d^)i] , (39) 
I 

which can be derived explicitly from Eq. (|36p . The modi- 
fied XY' model ([5^ possesses a second order phase tran- 
sition at the critical point^^ 

lim (/3) « 1.6. (40) 

The Kh 00 edge is trivial, being occupied by the su- 
perconducting phase with both the holon and the spinon- 
pair condensates present. 

The model alone was supposed^^ to possess an 
interesting link to the physics of correlated electrons be- 
ing able to describe certain exotic phases. The topolog- 
ical objects of this model are called "visons" which are 
fractionally charged excitations. In the language of the 
Abelian two-Higgs model the vison coincides with the 
spinon vortex, while the holon vortex turns into a XY 
vortex, ip becomes the field of the so-called "chargon" 
particle. The SG phase corresponds to the fractionalized 
phase where visons are absent and chargons are free par- 
ticles. In the SM phase the visons are condensed, and 
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chargons are confined. The SC phase corresponds to a 
superfluid state where both visons and XY vortices are 
dilute and chargons are free. 



B. The (/?, Ks)-faces 

The coupUng between the holon field and the gauge 
field is controlled by the parameter k^. Below we con- 
sider the influence of the spinous on the phase diagram 
considering the limits of large {k^ — > oo) and small cou- 
plings {kh 0). 

1. The Kh = face: Q — 2 compact Abelian Higgs model 

On the Kh — Q face the holon condensate vanishes, 
& = 0, and the spinon-pair condensate A is described by 
the cAHMq=2 model: 

S'cAhmq^s = -/3^ cos 6'p - Ks^cos.{d^ps + 26*); .(41) 
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FIG. 5: The phase diagram of the compact Abelian Q = 2 
Higgs model (|4ip on the k/i = face. 



The phase diagram of this model in the (/3, Ks)-plane is 
well known and the behavior of the spinon-pair conden- 
sate can be deduced from results of Ref. [27l[39H40l Before 
doing so, let us consider the edges of the phase diagram. 

The edges = and /3 = do not possess any phase 
transitions because they are described, respectively, by 
the cQED model ([5^ and by a trivial model. 

The edge (3 ^ oo corresponds to the XY model, since 
in this limit the constraint d0 = is imposed. The con- 
straint is resolved by setting = da, with subsequent 
identification 



shown by the dotted line in Figure [3 The SM phase is 
located at small /3 and/or small k^, while the SG phase 
is residing in the large /3/large Kg corner. The transition 
line corresponds to the onset of spinon-pair condensation 
A. In the whole diagram the holon condensate is absent, 
6 = 0. In addition also the limiting models {XY , cQED, 
gauge Ising and trivial cases) together with the critical 
points ^, Eq. ((Ml), and , Eq. ([381), are indicated. 



2. The Kh — > oo face: XY model 



+ 2q; — 27:1 = (z tt, tt] 

with / G Z. Then we obtain the XY action ([55)1 by set- 
ting K^^ = Kg. The XY model describes the superfluid 
behavior of the spinon-pair condensate with the XY crit- 
ical point ([33]), K^^. 

The edge Kg oo corresponds to the gauge Ising 
model ((37)) because in this limit the constraint 

dips + 20 ~ 2-Km , m g Z , 

is automatically imposed. The constraint is resolved as 

6 = nm — d(y5s/2 + 2T:n , n S Z , 

leading subsequently to 

d9 ~ Trdm , cos AQ = cos Trdm = up , 

where cr; — (—1)™' is the Z2 gauge field attached to the 
link I. Then Eq. (gT]) reduces to Eq. and the edge 

Kg ^ 00 possesses a critical point at /3c = /3f^ which 
marks the second order phase transition of Ising type. 

The critical points at the — > 00 and Ks — > 00 edges 
are connected by a second order phase transition line as 



On the K/i — > cxo face the model ([8|) reduces to the XY 
model ([55)) with a hopping parameter k = Kg and — 
—2iph -(- ^fg. The model controls the superfluid behavior 
of the spinon-pair condensate A. Due to the constraint 

9 ^difh + 2TTl, / e Z , 

the holon vortices are suppressed and therefore 6 in 
the whole (/3,Ks)-plane. The phase diagram - shown in 
Figure [6] - is divided by a second order XF-like transi- 
tion hue parallel to the (3 axis at Kg^dP) = shown as 
dotted line. This line separates the FL phase (with con- 
densed spinon vortices and A = 0) at /ts < from the 
SC phase (with suppressed spinon vortices and A ^ 0) 
at Kg > K^^ . In the whole diagram the holon field is 
condensed, 6^0. 



C. The (Kh, Ks)-faces 

Now we consider the effect of choosing the extreme 
limits of very strong and very weak gauge couplings, /3 = 
and /? — > 00, respectively. 
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FIG. 6: The phase diagram of the XY model (|33p together 
with the limiting models at the edges on the Hh ^ oo face. 



1. The (3 — face: an uHralocal two-Higgs system 

On the /3 = face we obtain a two-Higgs system in- 
teracting ultra-locally via a non-propagating gauge field. 
The one-link action of the model is given by 



/ D6'exp| ^ KkCOs{Aipk + qkO)i^ ■ (42) 



line) which starts at Kh^c = , Eq. (|40l) . at the Kg 
OO edge and ends at Ks.c = ^-^^ , Eq. ([5^ . for — ^ oo. 
At these two edges the two-Higgs model (El) is reduced to 
a modified ([5^ and a usual ([55]) model, respectively. 
The SM-FL phase appears actually as the FL phase (6 ^ 
0) at large Kh, and the SM phase (6 0) is realized at 
large Kg, (the structure of the SM-FL phase is plotted in 
Figure 121). 




Koo 



FIG. 7: The phase diagram of the ultralocal two-Higgs sys- 
tem on the /3 = face. 



k—h.s 



In order to get this part of the weight for the limiting 
model we put /3 = in Eq. ([8]) and keep the integration 
over the gauge field 9. We expand the two factors of the 
exponentiated one-link action in Eq. (j42p in a Fourier 
series, 

for g/i = 1 and Qs = 2. Here /„(k) is the modified Bessel 
function of 7i-th order. Substituting Eq. in Eq. (gH) 
and performing the integration over 9 we obtain the con- 
straint Uh + 2ns — 0. Setting Hs — — n/i/2 = n, we get 
(up to an inessential factor in front of the sum) 

where we used the property /„ = Note that the 

quantities n, diph, dips in Eq. (|44| are defined at the 
same link I. 

In the (K/i,Ks)-plane one has two phases: the SC phase 
with non-zero condensates A and b in the large k/i /large 
Ks corner and a SM-FL phase in the remaining part of the 
phase diagram. The phase diagram is shown in Figure [T] 
The SC and SM-FL phases are separated by a second 
order XY-type phase transition (indicated by a dotted 



2. The l3 —> OO face: two decoupled XY models 

Finally, on the /? ^ oo face the system 1^ reduces 
to two decoupled XY models describing the holon and 
spinon-pair superfluids. This fact is readily seen from 
the partition function ([8]) : large (3 imposes the constraint 
d9 = which is resolved, as usual, hy 9 — da + 2?™, n £ 
Z. Substituting this solution back to ^ and performing 
the redefinitions of the phases, iph + ct ^ y^h and tp^ + 
2a — > ifis, one gets 

S2XY = -Kh ^ cos{diph)i - i^s^ coa{dips)i ■ (45) 

The phase diagram in the (k/i, Ks)-planc includes all dis- 
cussed phases (SM, SG, SC and FL) as shown in Figure[51 
The phases are separated by the straight solid (conden- 
sation of the holon b) and dashed (condensation of the 
spinon-pair A) phase transition lines at Kh^c = 1^'^^ and 
Ks,c = K^^, Eq. respectively. 

D. The interior of the 3D phase diagram 

Knowing the limiting cases shown in Figures [3][5] allows 
us (following Ref. i2l} ) to reconstruct the interior of the 
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FIG. 8: The phase diagram on the /3 — > oo face of the two 
decoupled XY models (|45|) . 



three-dimensional phase diagram as shown schematically 
in Figure [SI In other words, the six faces of the 3D 
cube Figure ini correspond to the 2D diagrams plotted in 
Figures [313 




FIG. 9: The qualitative 3D phase diagram of the cA2HM 
(reproduced from Ref. [2ll). 

The shaded section in Figure [9] represents the phase 
diagram in the (k/j, Ks)-plane at fixed finite gauge cou- 
pling f3 in the region /?f^ < /3 < oo. A priori there are 
two possible views to expect of the internal section of the 
phase diagram. One of the options is plotted schemati- 
cally in Figure [TUl (a). The holon and spinon-pair con- 
densation lines are getting slightly curved with respect to 
the limiting {(3 oo) case shown in Figure [3 the solid 



(dashed) line, which marks the condensation of the holon 
b (spinon-pair A), gets generally shifted towards larger 
values of (ks)- The holon condensation line starts at 
large values of Kg at the point 



K/i,c(/3) > Kh,c{°°) 



This point is marked as the point C in Figures [3 and 
[TUl (a). With Kg becoming smaller the 6-condensation 
line meets the A-condensation line at the point F, and 
they continue together till the point G, at which these 
transition split again. Then the 6-condensation line con- 
tinues alone and eventually stops at an endpoint E in the 
interior of the diagram. A projection of this line to the 
edge Kg = (marked in Figure [TO] (a) by the dotted line 
E-A) might eventually be visible as a percolation tran- 
sition. The A-condensation transition, denoted by the 
dashed line D-F-G-B, does not have an endpoint. 
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FIG. 10: The 2D section of the 3D phase cube at a fixed 
gauge couphng p: (a) Pf < jS < oo and (b) < /3 < /Jf. 
The meaning of the points A,B,... is explained in the text. In 
both Figures in the middle of the diagram the transition lines 
of b- and A-condensation join along a finite segment. Figures 
(c) and (d) represent the respective alternative scenario with 
the transition lines intersecting instead of piecewise merging 
together. 

The two-dimensional section qualitatively changes as 
/3 gets smaller than the critical coupling /3f^ of the gauge 
Ising model. An example of such a section is shown 
in Figure [TUl (b). At sufficiently strong gauge coupling, 
(3 < /3f^, the points C, D and F merge together into 
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the point denoted now as H and the SG phase pocket 
disappears. The b- and A-condensation Unes begin to 
run together from the point H till the point F' (formerly 
point G in Figure [TO] (a)) in the interior of the phase di- 
agram, where the lines split again: the A-condcnsation 
goes to the point B', while the 6-condensation line runs 
towards the = edge but ends at a new endpoint E. 
As P decreases further, Figure [TO] (b) gradually evolves 
into Figure [T] at /? = the endpoint E becomes finally 
the point H. As a result, the 6-condensation transition 
disappears completely, as it is plotted in Figure [7] 

The alternative to the above scenario - which cannot 
be excluded by analytical means - is represented in Fig- 
ures [To] (c) and (d). The difference is in the mutual 
behavior of the b- and A-condensation lines in the in- 
terior of the phase diagram. In Figure [TO] (a) these tran- 
sitions piecewise join along a common segment F-G. In 
Figure [TO] (b) the common segment is H-F'. The alterna- 
tive of Figure [TO] (a) is plotted in Figure [TO] (c) : instead 
of merging along the F-G segment the transition lines in- 
tersect in the point F. The alternative of Figure [TO] (b) 
is plotted in Figure [TO] (d) : instead of having the com- 
mon segment H-F' (in Figure (b)) the b- and A-transition 
lines have only a common starting point H and separate 
immediately. 

To discriminate between these scenarios numerical sim- 
ulations have to be used. This is presented in the next 
section. 



IV. NUMERICAL RESULTS 

A. Observables 

In order to clarify the structure of the phase diagram 
we have performed a numerical study of various gauge- 
invariant quantities. One potentially sensitive quantity is 
the action of the model, or any (gauge-invariant) part of 
it. This fact is easy to understand since the action gov- 
erns the dynamics of the whole model. Using a familiar 
notation, the action of our lattice model is given as 

S = -13 Sp - Kh Sh - Ks Ss (46) 

with 

Sp — COs{9x,f_i + dx+fi,!^ ~ dx+o.n — f^x.v) , 

Sh^^ cos{Lph{x + fi)- fhix) + Ox,f,) , (47) 

S's = ^ COs{ips{x + fj.) - ips{x) + 2 9x,t_i) ■ 
x,^ 

Here x denotes the sites of the 3D lattice, nearest neigh- 
bors are separated by a lattice spacing a, and jx is the shift 
vector in the /i direction. The compact lattice gauge field 
angles 0^{x) G (— 7r,7r] live on the links (bonds) between 
sites X and x + p,, the holon fields (ph{x) (phase angles of 



a singly charged scalar field) and the spinon-pair fields 
ips {x) (phase angles of a doubly charged scalar field) are 
defined on the sites. The coupling (3 — l/g^ is the inverse 
gauge coupling squared, and Kh and Kg express the cou- 
pling between Higgs and gauge fields, respectively. In the 
following Kh and called "hopping parameters" . To 

study this model by means of Monte Carlo we use stan- 
dard Metropolis updates for all three kinds of fields. 

In order to characterize the different phases of the 
model we consider the following "thermodynamical" ex- 
pectation values (related to the derivatives of the loga- 
rithm of the partition function with respect to the cou- 
plings): 

{EP) = {^^Sp) 

{En) = {^Sh). (48) 

{Es) = (^^.> , 

called plaquette, holon link and spinon link expecta- 
tion values. Here Np/i^/g is the number of plaque- 
ttes/links/sites on the finite lattice {Np = Nl in three 
dimensions), (• • • ) denotes the ensemble average over 
configurations. In addition, the susceptibilities of these 
quantities (related to the second derivatives of the loga- 
rithm of the partition function with respect to the cou- 
plings) have been considered: 

X., = Np({El) - {Ep)^) , (49) 
X.,,, = Nl({EIJ - {Eu,s)^) . 

The simplest characteristics of a topological defect is 
its density. Using the notations of Section [IT] the mono- 
pole and the vortex densities are defined as, respectively, 

Praon = ^^\* j\ , PKs ^ ^^\* <JhA ■ (50) 

*C3 ^ *C2 

Here *C3 are the sites of the dual lattice dual to the cubes 
of the original lattice, and *C2 are the links of the dual 
lattice dual to the plaquettes of the original lattice. The 
monopole charge is defined in the standard way, 

J = ^d[d0]2^ , (51) 

where [■ • ■]2^/(27r) denotes the integer part modulo 27r. 
Note that the plaquette angle dO lies in the range — tt to 
TT plus or minus integer multiples of 27r. Thus, effectively 
using the Gauss's law, one considers a forward cube from 
a lattice position x built up by six plaquettes, defines the 
integer-valued so-called oriented Dirac strings n passing 
through these plaquettes and sums over these integers as- 
signed to the outward pointing oriented strings. This sum 
defines the monopole number in the cube corresponding 



to a point of the dual lattice. Following Ref.|4l|, the holon 
and spinon vortex currents are defined as 

(Tfe = -^(d[d(^ + gfc(?]2^ -<7fe[d0]2^), k^h,s. (52) 
Ztt 

These integer- valued oriented currents pierce a given pla- 
quette corresponding to a link of the dual lattice. To- 
gether with pnion and ph.s we measure also the corre- 
sponding susceptibilities 



(53) 



In this Section we report on Monte Carlo studies of the 
phase transitions within two-dimensional cross-sections 
of the whole phase diagram characterized by certain fixed 
values of the (inverse) gauge coupling /3. These planes are 
parameterized by the two hopping parameters Hh and Ks- 

It is important to note that we have to distinguish be- 
tween the cross-sections lying above and below the criti- 
cal value /3f^ « 0.7613, corresponding to the phase transi- 
tion of the pure gauge Ising model. As we have discussed 
above, the qualitative structure of the phase diagram is 
different above and below this phase transition, obtained 
by moving the shaded square in Figure [9] up and down. 
The SG phase (the phase with zero holon condensate and 
non-zero spinon-pair condensate) is absent in the f3 < (3^ 
(strong gauge coupling) region contrary to the (5 > (3^ 
(weak gauge coupling) part of the phase diagram. 

In our previous work^ we restricted ourselves to the 
larger (3 region and presented results mainly for (3 — 1.0. 
In order to obtain the gross features of the different 
phases, we perform the initial studies on 16^ lattices, 
focussing at different fixed values of (3. Based on the ex- 
pectations as described in the previous Section IIIIl we 
have used a dense grid of points spanning the (k/i,Ks) 
hopping parameter plane over an range where we expect 
a nontrivial behavior of the model. 



B. Phase structure at weak gauge coupling 

According to our analysis presented above, in the weak 
coupling region, (3 > /3f^, the phase diagram contains all 
four phases. The summary of our results on the phase 
structure of the model is presented in Figures [TT]fT^ 
These Figures show the two-dimensional cross-sections 
- explored, respectively, at /? = 2.0, /3 = 1.5 and 
f3 = 1.0 values of the gauge coupling - of the full three- 
dimensional phase diagram. 

The non-trivial signals in the measured susceptibilities 
indicate the existence of transitions between all four dif- 
ferent phases known already from our discussion of the 
weak coupling limit /3 — > oo, presented in Section IIIIl 
One can clearly observe that at the largest measured (3 
the four different phases meet (within our numerical res- 
olution) in a single crossing point. The numerically ob- 
served picture is very similar to the expected behavior 




FIG. 11: The two-dimensional cross-section of the three- 
dimensional phase diagram (Figure [9]) at /3 = 2.0, extracted 
from a 16^ lattice iisinff snscentibilities of thermodvnamical 
obs 
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FIG. 12: The same as Figure [TTl but for P = 1.5. 
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FIG. 13: The same as Figure [TT] but for /3 = 1.0. In addition, 
the four phases and positions of the points A, B, C and D of 
Figure |9] are explicitly indicated. 
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in the limit /3 — > cx) as shown in Figure [8l The hori- 
zontal (blue) and vertical (red) phase transitions are of 
second order, and both belong to the XY universality 
class. Note that here we do not discriminate between 
the ordinary and "inverted"— universality classes of the 
XY transitions as they are related by a turnover of the 
coupling axis, k —^ 1/k, in the course of a duality trans- 
formation. 

With decreasing /3, the gauge coupling becomes 
stronger and the phase picture changes (see Fig. [T3|) . 
Consider first the horizontal transition line which marks 
the condensation of spinous. At vanishing holon cou- 
pling, Kh = 0, the phases are separated at a certain 
"lower border" value (point D) of the critical spinon cou- 
pling, k"'*(k/i 0,/3), which turns out to be a rapid 
function of the gauge coupling /3. Indeed, as the cou- 
pling (3 decreases, the critical spinon coupling k"'* sub- 
stantially increases. On the other hand, at large values of 
the holon coupling k/j ^ oo, the "upper border" (point 
B) of the critical spinon coupling K"''(Kft, oo,/3) is 
practically insensitive to the variation of (3. 

The features of evolution of the vertical transition line 
- which corresponds to the condensation of the holon 
pairs ~ is similar to the evolution of the horizontal line. 
The lower endpoint of the line (point A) evolves moder- 
ately while the upper endpoint (point C) practically does 
not move at all. As we have already pointed out in Ref 
[2ll . the indicated line for small at P = 1.0 near the 
point A does no longer belong to a real phase transition, 
but characterizes a percolation transition from the SM 
to the SL phase. This is in agreement with our findings 
at lower /3's discussed below. 

It is very interesting to find, what happens with these 
two transition lines in the interior of the phase diagram. 
As it is already clear from Figure [121 the two transition 
lines do not simply cross each other in an isolated point. 
As the inverse gauge coupling (3 gets lower, in the middle 
of the phase diagram the transition lines become closer 
and closer to each other in a particular interval of the 
(ks,k/i) coupling space. The further decrease of the cou- 
pling (3 leads to a qualitative change of the picture in 
the interior, as it is indicated in Figure [T31 a,t (3 = 1.0 
the transition lines piecewise join into a single line in a 
certain {Kh,K,s) region. In this region we observe non- 
trivial signals for a first order transition in all measured 
quantities to be discussed in the next subsection. 

In addition to the "thermodynamical observables" 
the "topological quantities" (densities of defects) 
also show signals of a phase transition. The density of 
monopoles Pmon (at [3 = 1.0 on a 16"^ lattice) is plotted in 
the upper panel of Fig. [Til over the (K;i,Ks)-plane. With 
increasing hopping parameters Kh or Kg the monopole 
density gets suppressed. As shown in the lower panels 
of Fig. [TH the density of the holon vortices pt (spinon 
vortices ps) significantly drops down with increasing 
(or Kg, respectively). This behavior is not unexpected 
because, as one can see from the vortex action ([?7l) - ([5ni) . 
the larger the hopping parameter, the bigger the vortex 




FIG. 14: The monopole density (top) and the vortex densities 
(bottom) at /3 = 1.0 and 16^. The lower panel: the density of 
holon vortices (left) and of spinon vortices (right). 



mass. Therefore, the increase of a particular hopping pa- 
rameter must suppress the density of the corresponding 
vortex. As for the monopoles, the increase of either of 
the hopping parameters should suppress the monopole 
density because the monopoles are connected by vortices 
according to Eq. ([20l) (see the example in Figure [J). The 
increase of tension (mass) of either of the vortices leads 
to the confinement of the monopoles into magnetically 
neutral monopole-anti-monopoles states, and, as a result, 
this leads to the suppression of the monopole density as 
we see in Figure [T4l ftop). The maximum in the monopole 
density is seen where both vortex densities are non-zero. 



C. Strengthening of the phase transition: 
"2" -\- "2" = "1" 

The structure of the phase diagram at /3 > /9f' gives us 
the unique possibility to study the phenomenon of merg- 
ing (to be distinguished from crossing!) of two different 
phase transitions in a finite region of the coupling space. 

In order to clarify the nature of the phase transitions, 
we studied the volume dependence of the average of the 
plaquette and of the both link terms in the action ^ 
as well as their respective susceptibilities in different re- 
gions of the phase diagram. Figure [151 (top) shows the 
jump developing in the holon link vs. at fixed Ks 
with increasing volume. This jump is a clear signal of 
the first order nature of the transition (it is actually ob- 
served in all parts of the action!) in the crossing region 
of the transition lines. Note, however, that a lattice with 
a size 32"^ turned out to be too large to tunnel for the 
selected values, even within 5 x 10^ Monte Carlo it- 
erations. The reason is that the free energy barrier for 
such a large lattice is too high. In Figure [TSl (bottom) 
we present typical two-state signals (here only for a 12^ 
lattice) of all three terms of the action at {Kh,K,s) close 
to the transition. The signals are strong along the di- 
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FIG. 15: Signatures for the first order transition at /9 = 1.0: 
the holon link expectation value vs. Kh for different volumes 
at Ks = 1.4 (top), and (unnormalized) histograms of the parts 
of the action in the pseudocritical region on a 12'' lattice (bot- 
tom). 



ticipated for the limit Ks —>^ 0. 

Due to the weakness of the signals we were not able 
to check, whether for larger but finite values of /3 (corre- 
sponding to weaker gauge coupling) there is still a line of 
first order transition somewhere in the region of the phase 
diagram close to the crossing point. In fact, as (3 increases 
the region in which the two transition lines join, tends to 
shrink. Therefore the properties of the bulk observables - 
probed by variations of the Kh hopping parameter - pro- 
vide a weak signal if the sequence of measurements does 
not pass through the right (crossing) point at the right 
(corresponding to a maximal variation of the bulk ob- 
servables) angle in the coupling space. Anyway, the first 
order phase transition should inevitably disappear from 
the crossing region in the region (3 —>■ oo as we surely 
know from Figure [51 

Increasing further the strength of the gauge coupling 
(decreasing /3 below (] = 1.0) one should finally approach 
the critical value of the gauge Ising model, /3f^ « 0.7613, 
below which the structure of the two-dimensional phase 
diagram qualitatively changes. In order to see what hap- 
pens there we studied in detail the cross-section of the 
phase diagram at /3 = 0.8 which is already very close to 
the critical value. In Figure [TBI we present the plaquette 
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red transition line between the SM and SC phases. The 
two-state signals of the volume-averaged plaquette and 
volume-averaged holon link term become very weak when 
one goes to smaller along the horizontal dark-dotted 
(blue) line. Therefore in the crossing region between the 
two phases the strength of the phase transition is en- 
hanced compared to the strength of the individual (hori- 
zontal and vertical) transition lines. We discuss the order 
of these transition lines below. 

It is known that for k/i — > the transition vs. Hs is 
of second order. Similarly, for large the transition is 
most likely also of second order, again in the XY uni- 
versality class. We found that at Kh = 2.0 already for 
the largest volumes 40^ — 48^ of our study the increase 
of the spinon link susceptibility stops as function of the 
lattice size. This is the behavior expected for the XY 
model at — > oo^. Concentrating on two Kg values 
outside the crossing regime, where one of the transition 
lines (the light-dotted [red] one) runs vertical, we observe 
that there is no thermodynamic transition vs. Kh for the 
smaller Ks. This is in agreement with what could be an- 



FIG. 16: Plaquette expectation value in the {Kh, Ks)-plane at 
P = 0.8 on a 16^ lattice. 

expectation value (Ep) in the {Kh, Kg) -plane. The white 
areas in the Figure correspond to uninteresting regions 
in the plane of the couplings at simultaneously large Kh 
and Ks values. From this Figure we may expect that the 
phase transition line corresponding to the spinon conden- 
sation is parallel to the Kh axis for not too small Kg . The 
same effect can also be seen in the landscape of the av- 
erage spinon action (Es) (not shown). Analogously, the 
holon condensation line is parallel to the for not 

too small Kh as can also be seen in the landscape of the 
average holon action (Eh) (not shown). 

The maxima of the susceptibilities of the different ac- 
tion contributions demonstrate in more detail the exact 
location of the phase transition curves. To obtain accu- 
rate pseudocritical couplings we performed several high 
statistics runs near chosen points indicated in Figures [T71 
changing one of the hopping parameters to cross the ex- 
pected phase transition line. 
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All obtained Monte-Carlo histories of measurements 
have been evaluated together using a multihistogram 
reweighting procedure^. The combined histograms of 
observables at a chosen pair of couplings , Ks (interpo- 
lating among the K;t,Ks-grid of the used data sets) are 
obtained by reweighting from the closest data points in 
the grid. They contain all information necessary to pre- 
cisely locate the phase transition. In particular, using 
such a combined histogram interpolating in the range of 
the chosen hopping parameters we were able to identify 
the pseudocritical value where the reweighted histogram 
exhibits a maximum in the susceptibility of the corre- 
sponding variable. 

To estimate the error in determining the critical cou- 
plings, we blocked our original data from several cou- 
plings into blocks similar to a jackknifc method and con- 
structed individual "subhistograms" at the chosen cou- 
pling pair for those data subsets. Using those subhis- 
tograms to find the maximum of the susceptibility, dif- 
ferent critical couplings have been found which allowed 
to estimate the accuracy of the location of the critical 
coupling using all available data. 

The results of this detailed investigation is shown in 
Figures 1171 We used thermodynamical variables in the 



the susceptibility of the spinon vortex density. The ob- 
served behavior can be interpreted as indication for the 
expected second order phase transition. 




FIG. 17: The phase transition lines in the (Kh,Ks)-plane ob- 
tained by searching for peaks of various susceptibilities using 
multihistogram reweighting for (3 = 0.8 on the 16^ lattice. 

reweighting as indicated in the legend. Most errors in 
the pseudocritical couplings are smaller than the symbol 
sizes, the lines are drawn to guide the eyes. 

The individual observables show common signals in 
certain parts of the phase plane. The region of the com- 
mon transition line has been increased compared to re- 
sults at larger /3. In addition to the gross structure the 
susceptibility of the plaquette action shows a weak but 
measurable signal for the continuation of the "horizon- 
tal" line at large k™* « 8.0 for very small Kh- 

This has been studied in more detail at small Kh = 0.1 
in Figure [181 Weak but relatively sharp signals are ob- 
served in some of the susceptibilities shown as functions 
of the spinon hopping parameter . Note that contrary 
to the spinon link susceptibility, a signal is present in 
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FIG. 18; The susceptibilities of thermodynamical and topo- 



logical quantities at /3 

Kh = 0.1. 



0.8, and 16 as function of Ks at 



The emerging picture is in agreement with our expec- 
tation, that the critical value of the complementary Kg 
should take an infinite value as the value of /3 approaches 
the critical coupling /3f^ of the gauge Ising model. In- 
deed, there is no transition at finite Ks for the (3 = 2 
Abelian Higgs model at /3f^. 

Figure [TT] shows the location of the phase transition 
lines in a large region of the (K^,Ks)-plane. In particular 
we are convinced that the holon and the spinon-pair con- 
densation lines join for strong gauge couplings just above 
the gauge Ising model's critical value. 

In order to demonstrate that at /? = 0.8 the phase tran- 
sition in the merging region is of first order we have con- 
sidered three different fixed Kg values. Kg = 2.0, 4.0 and 
6.0. In Figure [19] we present for Kg = 4.0 near criticality 
nice two-state signals visible in various observables, indi- 
cating that the systems jumps from one metastable state 
to another. This is a clear characteristic of a first order 
phase transition. The signal is strongest near the center 
value Kg = 4.0, and then becomes weaker for smaller and 
larger Kg values (not shown) where - according to our 
expectations - the nature of the transition gets closer to 
a second order transition. 

The tunnclings between different vacua in the region 
of the first order phase transition are clearly visible for 
small lattice volumes (i.e., 16"^). They can also be seen 
in the (unnormalized) distribution of the holon energy 
term shown in Figure [^Dj for different lattice volumes and 
various values of the holon hopping parameter. 

The typical features of the first order phase transition 
can also be observed, for example, in the behavior of the 
holon link expectation value. This observable - shown in 
Figures [21] as a function of Kh for the same three chosen 
values of the hopping parameter Kg - clearly develops a 
jump at the phase transition point as the volume of the 
lattice increases. Thus the model has a non-vanishing 
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FIG. 19: Two-state signals in the Monte Carlo history of 
several observables near the center (ks = 4.0) of the merging 
line of the two transition line as observed at /3 = 0.8 on a 16^ 
lattice. 




FIG. 20: The two-state signal in the histogram of the holon 
link expectation value Eh = Sk/Nl at /3 = 0.8 and Hs ~ 4.0 
with varying lattice volume reweighted to the transition value 
Kh at maximal susceptibility. 



latent heat at the phase transition point characterizing 
a first order phase transition. The signal of the first or- 
der is most clearly visible at Kg = 4.0, while it becomes 
weaker at larger and smaller values of Kg in agreement 
with our observations made with the help of the Monte 
Carlo histories. 

In order to confirm the absence of the phase transition 
as we expect from either Figure [TU] (a) or Figure [TU] (c) , 
we investigated the model numerically for a much smaller 
value of the spinon hopping parameter, Kg = 0.1. The 
crossover nature of the transition is clear from the sus- 
ceptibilities of the thermodynamical and topological ob- 
servables (in Figure : we observe either no suscepti- 
bility signal or broad maxima at different values of k^. 
Of course, a remnant of a vortex percolation transition 
- expected at small Kg - cannot be ruled out by our in- 
vestigations and has to be studied using cluster analysis 
techniques. 

Thus, we have strong reasons to conclude that the 
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FIG. 21: Holon link energy expectation value as a function 
Kh for various volumes at /? = 0.8 for Ks = 2.0 (top), Ks = 4.0 
(middle) and Ks — 6.0 (bottom). 



phase transition scenario at relatively strong gauge cou- 
pling, but still above the critical /3f' of the gauge Ising 
model, resembles more one of the anticipated scenarios, 
the phase structure of Figure [TO] (a) rather than the clos- 
est alternative scenario plotted in Figure [TO] (c). 

Summarizing, we have clearly observed the strength- 
ening of the phase transition in the merging region. The 
strengthening happens due to activity of the internal 
gauge field in the weak coupling region. 
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FIG. 22; The susceptibilities of thermodynamical and topo- 
logical quantities at /3 = 0.8, and 16"^ as function of kh at 
Ks = 0.1. 



D. The phase structure at strong coupling 
below I3f 



In this Section we identify the phase structure of the 
model choosing a really strong gauge coupling, /3 = 0.5 < 
/3f^. As in the previous subsections we first make an 
exploratory study of the phase structure for a relatively 
small lattice size, 16^. 

In Figures [23] we show the plaquette expectation value 
and the average densities of all topological defects in the 
(k^, Ks)-plane (for the topological densities the viewpoint 
is rotated for convenience). Using the gauge degrees of 
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FIG. 23: Expectation value of the gauge plaquette and the 
densities of the topological defects (monopoles, holon vortices 
and spinon vortices) at /3 = 0.5 on a 16"^ lattice. 

freedom, the approximate phase structure is most easily 
detected. 

It is interesting to discuss the density of the topological 
defects since they may provide us with the hint how the 
condensates behave in the different regions of the phase 
diagram (a condensate is suppressed if the density of the 
corresponding vortex is non-zero). The density of the 
holon vortices behaves qualitatively similarly to the one 



in the weak coupling case. Figure fT4l In contrast to this, 
the density of the spinon vortices is qualitatively differ- 
ent at strong and weak gauge couplings as a function of 
Kg and K^: while in the weak coupling regime the spi- 
non vortex density rapidly vanishes going to the large-K^ 
limit, the same is no longer true in the strong coupling 
regime. The monopole density is also different in both 
regimes: in the weak coupling regime the monopole den- 
sity is visibly non-zero only in the vicinity of the Kh , Ks 
origin, while in the strong coupling regime the monopole 
density remains almost unaffected towards large values 
of Kg. Note that similarly to Figure [HI the rule remains 
true that a non-zero monopole density is observed only 
in the region of the phase diagram where both vortex 
densities are non- vanishing. 

Figures[23]show approximately a structure of the phase 
diagram qualitatively consistent with our expectations 
plotted in Figure [TUI (b). 

In order to study the location of the phase transition 
line more precisely, we have studied the thermodynamical 
and topological susceptibilities interpolated by reweight- 
ing techniques as for /3 = 0.8. The result of a search 
for maxima of those susceptibilities is represented in Fig- 
ures [21] 
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FIG. 24: The phase transition line in the (Kfi,Ks)-plane at 
fixed /3 — 0.5 on the lattice 16'^ to be located by finding the 
susceptibility peaks of thermodynamical (top) and topological 
(bottom) variables. 



We note that the different variables mark a single line 
in the hopping parameter plane (within our resolution). 
Signals from the thermodynamical susceptibilities (upper 
Figure IM)) are found in the whole considered {nh, i^s)- 
range. In contrary, using the topological observables 
(lower Figure [24|) , only the maximum of the spinon vertex 
susceptibility signals the transition along the whole 
line. The monopole and holon vertex signals are seen in 
the common "vertical" line only. 

The volume dependence of the susceptibilities of the 
holon link (top) and the spinon link (bottom) along lines 
crossing the phase transition line horizontally at Ks — 
4.0 and vertically at Kf^ = 3.0 definitely rules out the 
possibility of a first order transition (Figure [^5]) . 
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FIG. 25: The holon link susceptibility vs. at fixed Ks = 4.0 
(top) and the spinon link susceptibility vs. Ks at fixed kh ~ 
3.0 (bottom) for various lattice volumes at = 0.5 together 
with the curves obtained by multihistogram reweighting. 

Similarly to the weak coupling case {3 > /3f ^ , a branch 
of phase transition is absent ranging to very small Ks = 
0.1, as it can also be seen from the behavior of various sus- 
ceptibilities, in Figure [26l Indeed, these susceptibilities 
develop non-sharp maxima at significantly different k^. 
Note that the existence of a remnant from a vortex per- 
colation transition cannot be determined by our calcula- 
tions. It has to be sought for by cluster analyses. 

Summarizing this subsection we mention in short, that 
we observe one main phase transition line which is most 
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FIG. 26: Susceptibilities of thermodynamical and topological 
quantities vs. kh at (3 = 0.5, Ks = 0.1 and 16^. 



presumably of second order. In addition, we see a very 
weak signal along the vertical line k/j ~ 1 for small Kg < 
1.0 intersecting the Ks = axis. However, this mentioned 
signal is not a thermodynamical transition. 

The observed picture at /3 — 0.5 is consistent with 
the phase diagram proposed in Figure [TU] (b). We 
clearly observe the common line H-B' indicating the A- 
condensation as signaled by the susceptibility maximum 
of the spinon vortex density and the behavior of the den- 
sity itself. The "vertical" section H-F' of the Hue H-B' - 
showing in addition the 6-condensation ~ is also present 
as indicated by the behavior of the holon vortex density. 

There is no hint for thcrmodynamically separated tran- 
sitions for the condensations. If a line E-F' with a point 
E exists, it could be related to percolation properties as 
can be conjectured from Figures and 



V. THE STRENGTHENING OF PHASE 
TRANSITIONS AND A POSSIBLE 
REALIZATION IN QCD 

This paper is devoted to the investigation of non-trivial 
properties of a model containing two Higgs fields coupled 
to a single compact gauge field. As most remarkable ob- 
servation we find the possibility that two weak (second 
order) phase transitions - marking the vanishing/non- 
vanishing of two different condensates - tend to join and 
constitute a stronger (first order) transition at sufficiently 
strong (but not too strong) gauge coupling. Although the 
model is related to strongly correlated electron systems 
in (2-fl) dimensions (at zero temperature), we conjec- 
ture that the observed effect is much more general such 
that it can also be realized in the field theory of strong in- 
teractions, Quantum Chromodynamics (QCD). The two- 
dimensional phase diagram that attracts the attention of 
lattice theorists and phenomenologists today, is spanned 
by temperature and non-zero baryonic density. 

The schematic view of the QCD phase diagram in 
terms of the baryon chemical potential fi and the temper- 
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ature T is shown in Figure [27l For a detailed review we 
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FIG. 27: Schematic view of the phase diagram of QCD with a 
suggested enhancement to a first order phase transition due to 
the mediating, strongly coupled gauge field. The explanation 
is given in the text. 



refer to Refs. I45ll46l . The diagram corresponds to QCD 
with two light quarks {u and d) and one heavier quark (s). 
The order parameters which are relevant in this phase di- 
agram are the chiral condensate, s-nd the diquark 
condensate, ('0V') • It is worth mentioning that due to the 
presence of three quark flavors the phase structure of the 
realistic QCD is richer compared to FigurelTTl The struc- 
ture of the superconducting phases is finer because there 
are three (instead of one) diquark condensates related 
to different flavors, and the phases with different combi- 
nation of non-zero diquark condensates may emerge. A 
thorough treatment of the QCD phase diagram in an ef- 
fective model of the quark matter, known as the Nambu- 
Jona-Lasinio (NJL) model^, can be found in Ref. '48'. In 
view of the qualitative character of our considerations we 
ignore the color, flavor and Dirac indices in the notation, 
and we do not discriminate between different quark fla- 
vors taking the schematic picture shown in Figure [27] as 
a reference point. 

At low temperature and density, strongly interacting 
matter is in the hadronic phase which is characterized 
by broken chiral and unbroken color symmetries, so that 
the quark condensate is non- vanishing while the diquark 
condensate vanishes. As temperature increases, matter 
goes over into the quark-gluon plasma (QGP) phase in 
which both condensates vanish. At small chemical po- 
tential the transition between these phases is probably a 
smooth crossover**^ rather than a real phase transition, so 
that in this region the difference between the two phases 
is marginal (and the vanishing of {tp^j) is not complete). 
However, the chiral condensate is drastically suppressed 
in the high temperature regime, what we formally ex- 
press in Figure [57] by (V'V') = 0. We adopt the same 
convention for other exponentially suppressed but non- 
vanishing condensates. The phases on both sides of a 
crossover are often said to be connected via an analyti- 
cal continuation. Consequently, in a strict mathematical 



sense, the chiral condensate does not take zero value at 
any point of the quark-gluon plasma phase. 

As the chemical potential increases, the crossover be- 
tween the hadron and the QGP phases turns into a sec- 
ond order transition denoted by a dot in Figure WJ\ The 
second order phase transition is an endpoint of a first 
order phase transition line that separates the hadronic 
phase not only from the QGP phase but also from one 
of the color superconducting phases, the so called 2SC 
phase, in which two colors are locked with the global 
flavor symmetries through the formation of a non-zero 
diquark condensate, 7^ 0. Note, that in the 2SC 

phase the chiral condensate is zerc"*^. 

The 2SC phase is presumably also separated by a 
crossover from the QGP phase and by a first order phase 
transition from another color superconducting phase in 
which all three colors are locked with all three flavor de- 
grees of freedom (this phase is known as "color-flavor 
locked", or CFL, phase). The CFL phase is the only 
phase where both chiral and diquark condensates are si- 
multaneously non-vanishing^^ . In Ref. 50 it is argued 
that the transition between the QGP and the CFL phase 
must be of first order due to active role of the gauge 
fields. 

Needless to say that the system of strongly interacting 
gluons coupled to dynamical quarks in (3+1) dimensions 
is much more complicated than a simplified gauge model 
of correlated electrons in (2-f 1) dimensions. However, 
there are some qualitative similarities between these the- 
ories as far as the structure of the phase diagram is con- 
cerned comparing, say FigurefTUlfa) with Figure [?7l The 
compact Abelian two-Higgs model may be realized in one 
of four different phases: Fermi liquid, spin gap, strange 
metal and superconducting phases which are character- 
ized by pairs of vanishing or non-vanishing condensates 
(6- or holon. A- or spinon-pair condensates, see also Fig- 
ure [T]). The QCD phases can also be characterized by a 
pair of quark and diquark condensates {{ipip), (ipip)). Ac- 
cording to Figure [T51 an enhancement of the strength of 
the phase transition may take place (at sufficiently strong 
gauge coupling) where both condensates simultaneously 
change from vanishing to non- vanishing. The same pat- 
tern is realized in QCD between the CFL phase and the 
QGP phase. The corresponding segment of the phase 
transition is the result of joining two "transitions" , the 
crossover transition between the 2SC color superconduc- 
tor and the QGP on one hand and the first order phase 
transition between the 2SC phase and the CFL phase. 

Our experience with the cA2HM suggests that such 
a merging of two weaker transitions - provided they 
separate the phase characterized by simultaneously van- 
ishing condensates from the phase with all condensates 
non-vanishing - may lead to an enhancement of the or- 
der compared to the strength of the individual transi- 
tion lines. Therefore we expect that at sufficiently large 
baryon density the finite-temperature transition between 
the CFL phase and the QGP phases should be much 
stronger than between 2SC and CFL superconducting 
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phases taking place at lower temperatures. 

As we have discussed above, the key ingredient of 
the transition strengthening in the effective model of 
the strongly correlated electrons is a gauge-boson medi- 
ated interaction between the species of the matter fields. 
As a result, the transition in different channels/species 
merge and the merger corresponds to a much stronger 
phase transition. The strengthening effect cannot be seen 
in QCD phase diagram studies using the effective NJL 
model of quark matter (see, for example, Ref. |48| ) be- 
cause that model possesses global symmetries only and 
the gauge field mediation is obviously absent. There 
are, in fact, indications^ that in QCD the (thermal) 
fluctuations of the gauge fields strengthen the first or- 
der phase transition between the QGP and CFL phases, 
which nicely matches with our qualitative expectation. 

VI. CONCLUSIONS 

In this paper we have studied in detail the phase dia- 
gram of the (2-1- 1)D Abelian Higgs model with two Higgs 
fields and one compact gauge field. At weak gauge cou- 
pling (large /3) the phase diagram contains two transition 
lines which run through the whole plane in the hopping 
parameter (k?i, Ks)-plane. The transitions are associated 
with the onset of vortex percolation, and with the appear- 
ances of holon and spinon-pair condensates. The pattern 
of non-vanishing condensates in different regions of the 
phase diagram allows to identify the Fermi liquid, spin 
gap, superconductor and strange metallic phases. 

With decreasing (3 above the critical coupling of the 
gauge-Ising model the strange metallic and Fermi liquid 
phases become analytically connected at small Kg, even- 
tually still being separated by a percolation transition. 

At extremely strong coupling (with [3 below the gauge- 
Ising model transition), however, one of the segments 
of the transitions disappears together with the spin gap 
phase, while the difference between the strange metal and 



the Fermi liquid phases becomes almost invisible since 
these phases are separated by a smooth crossover in this 
limit. 

The most intriguing effect - which was first suggested 
in Ref. fH] - is a strong enhancement of the strength of 
the phase transition. This enhancement appears in the 
regime of moderately strong gauge coupling in which the 
two second order transitions - corresponding to spinon- 
pair and holon condensation lines - join. The enhance- 
ment happens along the common segment of the phase 
transition and is reflected in lowering the transition or- 
der: two second order transitions become a single first 
order transition between the two phases where the two 
condensates are both vanishing and both non- vanishing. 
That strengthening is a result of the activity of the inter- 
nal gauge field which couples the dynamics of the spinon- 
pair and holon condensates. 

We suggest that a similar enhancement effect may 
also be realized in Quantum Chromodynamics at non- 
zero temperature and at finite baryon density. At suffi- 
ciently large baryon density the finite-temperature tran- 
sition between the (3-flavor paired) color superconduct- 
ing phase and the quark-gluon plasma phases should be 
much stronger compared with the transition between 2- 
fiavor paired and 3-flavor paired superconducting phases. 
The suggested enhancement is emphasized in Figure [571 
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